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Abstract 

This work describes two algorithms for computing the 
angular rate and attitude in case of a gyro and a Star Tracker 
failure in the Microwave Anisotropy Probe (MAP) satellite, 
which was placed in the L2 parking point from where it 
collects data to determine the origin of the universe. The 
nature of the problem is described, two algorithms are 
suggested, an observability study is carried out and real 
MAP data are used to determine the merit of the algorithms. 
It is shown that one of the algorithms yields a good estimate 
of the rates but not of the attitude whereas the other 
algorithm yields a good estimate of the rate as well as two 
of the three attitude angles. The estimation of the third angle 
depends on the initial state estimate. There is a contradiction 
between this result and the outcome of the observability 
analysis. An explanation of this contradiction is given in the 
paper. Although this work treats a particular spacecraft, its 
conclusions are more general. 

INTRODUCTION 

The Microwave Anisotropy Probe (MAP) Satellite was 
launched at 15:46:46 EDT on June 30, 2001 aboard a Delta 
11-7425-10 launch vehicle. After three phasing loops MAP 
passed by a swingby point where it was injected by the 
moon’s gravity towards the L2 point. On October 1, 2001, 
following a three month journey, MAP arrived safely at its 
permanent observing station near the L2 Lagrange Point, 1 .5 
million km from Earth, a quasi-stable position a million 
miles from Earth in the direction opposite the Sun. 

MAP used two two-axis gyroscopes to measure its angular 
rate vector. One input axis of the first gyroscope was 
aligned along the body x-axis and the other input axis was 
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aligned along the body z-axis. The input axes of the 
second gyro were aligned along the body y and z-axes. 
Therefore, if one gyro failed then the rate about two body 
axes is still measured; namely, about the z-axis and one of 
the other two. MAP also had two Digital Sun Sensors that 
measured one vector; namely, the direction to the Sun. In 
addition, the spacecraft also carried two Star Trackers. 
However, when MAP performed the phasing loops it 
passed through the Van-Alien radiation belts, where the 
star trackers were inoperative. 

The problem handled in this paper is as follows. Suppose 
that indeed one gyro failed and the star trackers were 
inactive when MAP passed through the Van-Alien belts. 
Was it still possible to reconstruct the full angular rate and 
the orientation of MAP with the one operating gyro and a 
measurement of the Sun direction? It should be noted that 
although this question concerns a particular mission, its 
implication is more general. 

We addressed this problem by defining two possible 
filtering models and performing a nonlinear observability 
analysis. Then we performed simulations of the spacecraft 
motion and ran the suitable filters to verify the results of 
the analysis. In the next section we describe the angular 
rate at which the problem is examined. This section is 
followed by a description of the examined filters. In the 
section that follows we present the observability analysis 
that was carried out. In the next section we present 
simulation test results and then, in the last section, we 
present our conclusions. 

Body Angular Kinematics 

During the critical part of the trajectory; that is, when 
MAP passed through the radiation belts the spacecraft 
was kept at an inertial hold where the angular rate was 
nominally zero. In this mode the spacecraft attitude 
control system acted as a regulator where the gyros 
nominally read no input. In order to examine our 
suggested filters under more stringent conditions we 
chose to examine the filters under non-zero nominal 
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Fig. 1: The Geometry of the Angular Motion Components 
of the MAP Spacecraft. 


angular rate profile. As a baseline rate profile we chose 
the MAP science collection rate profile at L2. At the L2 
observation station the angular motion of MAP consists of 
two rotations as described in Fig. 1. One rotation, at 360 
degrees an hour, is about the anti-Sun line, and the other, 
at 0.464 RPM, is about the satellite z-axis. The body z- 
axis is tilted by 22.5 degrees with respect to the anti-Sun 
line. The body angular rate and attitude definitions are 
presented in Fig. 2 in terms of three Euler angles. In order 
to describe the spacecraft motion and the associated 
measurements we define four coordinate systems. The 
first system is the Sun coordinate system (designated by 
sun). The center of this system is at the origin of the 
spacecraft coordinates, its z-axis points at the Sun 
direction and its x-axis, which is turning, is assumed to 
point at some arbitrary direction in space, which we 
define as zero y . A rotation by the angle — y about the 
z-axis yields a first interim coordinate system, which we 
denote by 1. A second rotation by the angle 5 about the 
y-axis of this interim coordinate system yields the second 
interim coordinates, which we denote by 2. Finally, a 
rotation by the angle Tj about the z-axis of the latter 
system results in the body coordinates, which we denote 
by b. For sake of simplicity, only the body system z-axis 
is shown in Fig. 2. MAP rotates about two axes. One 
rotation, at the rate of — y , is about the Sun line, which 
we defined as the Z sun axis, and the other, at the rate of f] , 
is about the body z-axis. Using Fig. (2) we obtain 
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Fig. 2: Coordinate Systems Involved in the Description 
of the MAP Total Angular Velocity. 


© x = ysin 8 • cos rj 

o y = -y sin 8 • sin r| (1) 

G) z = -ycos8 + 

and the transformation matrix from the Sun to body 
coordinates is: 


D 


sun 

b 


cyc5cq + sysri 
-cyc5sr[ + sycrj 
cys5 


-sycbcT] + cysri 
sycSsri + cycri 
-sysS 


— sSCTj 
s5s r\ 
c5 


( 2 ) 


where s denotes the sine function and c the cosine. The 
corresponding quaternion can be extracted best when 
using the newly derived algorithm. 1 The angular rates and 
angles are as follows: 

y = 1 rev /hour = l-27r/3600 rad/sec 
and (3. a) 

y = yt + 0 rad = yt rad 

f| = 0.464 rpm = 0.464 -2n/ 60 rad /sec 
and (3.b) 

T]=f]t + 0 rad =f]t rad 
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5 = 22.5° =22.5- ni 180 rad 


(3.c) 

Upon examination of the MAP trajectory, the Sun direction 
and thus the Sun coordinates can be considered stationary 
in inertial space relative to the rotation of MAP with respect 
to the Sun coordinates. Consequently, we may consider the 
rate components of Eq. (1) as inertial rates resolved in body 
coordinates. 

Preliminary Analysis 

We wish to know the relationship between the sun sensor 
measurements and the angles y , 5 and T) that determine 
MAP’s attitude. To do that we turn to Fig. 3 which depicts 
the orientation of the sun sensor boresight in body 
coordinates. The sun sensor reads the tangent of the angles 
a and p . We wish, to express the readings of these 
tangents as functions of the angles y , 5 and T| defined in 
Fig. 2. To do that we define (Z sun ) b and (Z sun ) sun . The 
column matrix (z sun ) b contains the components of the 
vector Z sun when resolved in the body coordinates, and, 
similarly, the column matrix (Z sun ) sun contains the 

components of the vector Z sun when resolved in the Sun 
coordinates. We denote the components of the former by 
a v , a„ and a 7 . It is easy to see that the components of 

x y z 


(^sun )sun 

are 0, 0, and 1 . Thus 



(Z )I = Ta a a 1 

'^sun-'b |_ x y z J 

(4.a) 

and 


(Zsun)Ln=[0 0 1] 

(4-b) 

Obviously 


(Z sun )b = D b un (Z sun ) sun 

(5) 


From the last three equations and Eq. (2) we obtain 


a x 


-s5cr| 

'o' 

a y 

= 

x x sSsrj 

0 

_ a z_ 


c5 

1 


hence 

a x =-s5cr| 
a y = s5sq 
a 2 = c5 

therefore 

a 

tan P = — = - tan 5 • crj 

a z 

a y 

tan a = — = tan 5 • srj 


( 6 ) 


(7) 


( 8 ) 


Note that tana and tan (3 are independent of y. An 
inspection of Fig. 2 reveals that this is an expected result. 
Since tan a and tan (3 are measured quantities, the last 
two equations contain only two unknowns, therefore we 
can solve for the angles 8 and T] . For this reason even if 

the measurements contain unbiased errors we can still 
find these two angles by proper filtration. In order to fully 
know the spacecraft attitude we need to know y too. If 
the gyros mounted on the Sun coordinates z-axis 
introduce no bias, then a proper integration of their 
measured outputs will yield y and thus we can obtain the 
spacecraft attitude from the measurements obtained from 
these gyros and the sun sensor. If the gyro measurements 
contain a bias then we ask ourselves whether the system 
is observable; that is, whether the spacecraft angular 
dynamics (Euler) equation may add information with 
which we can still observe y and thus completely 
determine the attitude. To answer this question we have 
to design an estimator and examine the observability of 
the system used by the estimator. This is done next. 

THE FILTERS 

Filter I 

Dynamics Model 

In this estimator (filter) we convert the sun sensors 
measurements to a vector measurement, which is then 
connected to the quaternion of the spacecraft. The 
dynamics model of Filter I was developed in. 2 Using 


q4 

-43 

42 

43 

44 

-4i 

-q 2 

4l 

44 

,-qi 

-42 

-43 


define the system dynamics as 


(O 

q 


r 1 [(I« + h)x] 

0 

CO 

+ 

I _1 (T-h) 


0 

_q_ 


0 



( 10 ) 


where [gx] is the cross product matrix of a general 
vector g, accounts for the inaccuracies in the 

modeling of the spacecraft angular dynamics, and w q 
accounts for modeling errors of the quaternion dynamics. 

Measurement Equations 

The measurements are the tangent of the angles a and 
P , measured by the sun sensor and shown in Fig. 3, and 

the z and either the x or the y component of the body rate 
measured by one of the spacecraft gyros. The measured 
quantities are contaminated by noise so they are not the 
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reference values but rather, a m , P m , CD z m and either 
CO x m or C0 y m , depending on which gyro failed. 

Sun Sensor Measurement 

Using the two measured angles we compute a unit vector 
in the assumed sun direction as follows J 


b m 


1 

V tan + tan +1 


tan P m 
tana m 
1 


( 11 ) 



Fig. 3: Definition of the Sun Direction 

In reality tan p m and tan a m contain noise; however in the 

filter we approximate the b m measured vector as the true 
vector to the Sun and a simple additive zero mean white 
noise; that is, 

b m = b + v b (12) 

where b is the nominal unit vector in the true direction to 
the sun. The logic behind this simplification is that if in tests 
the filter yields satisfying results then the approximation is 
justified. Let r denote the measured Sun-vector expressed 
in the reference coordinate system (this vector is taken from 
the almanac). The relationship between the two vectors is 
expressed by 


b = Dr (13) 

where D is the direction cosine matrix that transforms 
vectors from the reference to the body coordinates. From 
the last two equations we obtain 

b m= Dr + v b 0 4 ) 

It is well known that D is the following function of the 
quaternion elements 4 


D = 


? 2 “> 2 

qf ~q2 -q§ + q4 

2 (qiq 2 ~q 3 q 4 ) 

2(qiq 3 + q 2 q 4 ) 


2 (qiq? + q 3 q 4 ) 

J2 , 5 /I _2 

-qi +q 2 -q 3 + q 4 
2(q 2 q 3 -qiq 4 ) 


2(qiq 3 -q 2 q 4 ) 

2 (q 2 q 3 + qiq 4 ) 

2 2 2 2 
-q, -q 2 + q 3 +q 4 


...(15) 


when the last expression is substituted into Eq. (14) we 
obtain 


K 


(qf -ql + 44)4 +2(q,q 2 +4344)^ +2(qiq 3 - 4244*3 
2(q,q 2 -4344)4 +(-qf +qf -q? +44*2 +2(q 2 q 3 +4i4 4 *3 
2(q t q 3 + 4244 *1 + 2(q 2 q 3 - qi44* 2 + (“4? “42+43 +44 *3 


(16) 


+ V b 


or 


qf r l + 4l42 r 2 +4l43 I 3 - 42 r l + 4l42 r 2 - 4244 r 3 “434 

qi42 r i - qf r 2 + qiq4*3 + qiq2 r i + qk + q2q3 r 3 - 43444 

qiq3 r i - qiq4 r 2 - q? r 3 + 42444 + q2q3 r 2 - qf r 3 + qi^n 


+q 3 q 4 r 2 +qiq 3 r 3 +q 4 r, +q 3 q 4 r 2 -q 2 q 4 r 3 
-q3 r 2 +q2q3 r 3 -434 4 q +q 4 r 2 +4i4 4 r 3 
+q 2 q 3 r 2 +4343+42444 _ 4 i 44 r 2 +4443 


+ v b 


which can be written as 


(17) 


b_ = 


Qi r i + q 2 r 2 + q 3 r 3 

q 2 q -qi r 2 + q 4 r 3 
q 3 r i ~q 4 r 2 -qi r 3 


~q 2 r i + qf 2 _ q 4 r 3 _ q 3 r i + q 4 f 2 + qi r 3 
m + q 2 r 2 + q 3 r 3 ~q 4 r i ~ q 3 r 2 + q 2 r 3 
q 4 q + q 3 r 2 - q 2 r 3 qi r i + q 2 r 2 + q 3 r 3 


Let 


q 4 ii+q 3 r 2 -q 2 r 3 
-q 3 q+q 4 r 2 + qi r 3 
q 2 r i-qi r 2 + q 4 r 3 


...(18) 


H b (r,q) = 


q,r, +q 2 r 2 +q 3 r 3 
q 2 r, -q,r 2 +q 4 r 3 
q3>l -Cl4r2 -Cll r 3 


-q 2 r i +qjr 2 -q 4 r 3 
qi r i +q2 r 2 +q3 r 3 
q 4 r ! +q 3 r 2 -q 2 r 3 


-q3 r i+q4 r 2+qi r 3 q4 r t+q3 r 2-q2 r 3 
-q 4 ri -q 3 r 2 +q 2 r 3 -q 3 r, +q 4 r 2 +q,r 3 

4l4 +42 r 2 +43 r 3 q2 r i - qi r 2 + q4 r 3 


(19) 


then Eq. ( 1 8) can be written as 


or 


b m =H b ( r ,q)q + v b 


bm=[0 


H b (r,q)] 



( 20 ) 

( 21 ) 


Eq. (21) is the measurement equation associated with the 
sun sensor measurement. 
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Gyro Measurement 

Let us assume that the gyro which measures rates along 
the x and z axes continues to operate. The jc and z axes rate 
measurements can be written as 


oo 

© 


x,m 

z,m 



+ v,. 


+ V„ 


( 22 ) 


where C0 X and co z are the true rates, and V C1)X and V ffiz 
are scalar zero-mean white-noise measurements. Define 


H 


CO - 


1 0 0 
0 0 1 


then from the last two equations we obtain 


(23) 


CO 


x,m 


CD 


= [H. 0 2I4 ] 


z,m 


CD 

q 


(24) 


We denote the measurement vector by G* m where, clearly, 
(ojj =[cD xm (D zm ]- Similarly we denote the noise 

vector as by where = [v^ . Obviously, if 

the other gyro is the one that is still operating, then the * 
component in Eqs. (22 - 24) has to be replaced by the y 
component. 


The Combined measurement 
Combining Eqs. (21) and (24) we obtain 


b m 


°3x3 

H b (r,q)' 

V 


'V 

_H W 

°2x4 . 

_qj 


_ v co _ 


(25) 


CD 

q 

b 


"r'[(Ico + h)x] 

0 

o" 

CD 

±Q 

0 

0 

q 

[bx] 

0 

0 

b 


,.n , „ (27) 



'r'(T-h)" 


’ w .* 

+ 

0 

+ 



Dr 


. w b. 


The state vector in Eq. (27) includes the Sun vector b 
even though b is a measurement. Indeed, in order to 
formally compute the rate and attitude one can 
numerically differentiate the measured b x and b y , and 

use them in the following equations (which are the 
components of Eq. (26.c)) 


b = co b - © b 

X y x x y 

b = © b - © b 

y z x w x z 


(28) 


together with the gyro-measured CD Z , to solve for the two 
unknowns CD X and CD y . Then use these CD X and CD y and 

the measured CD Z as ‘measurements’ in an estimator 
whose dynamics model consists of the first two matrix 
equations in Eq. (27). However, this method required 
numerical differentiation, which introduces noise. To 
avoid differentiation we apply the estimation approach . 5 
The latter can be accomplished by adding Eq. (26.c) to 
Eq. (10), the differential equations of the variables that we 
want to estimate, thus forming Eq. (27). At this point a 
Kalman filter can be utilized to estimate those variables. 


This is the combined measurement model for simultaneous 
measurement of the operating two-axes gyro and the sun 
sensor. 


Filter II 

Dynamics Model 

In this filter we make use of the kinematics rule governing 
the rate of change of the Sun vector. It is known that 

Dr = b + co x b (26.a) 


hence 


b = [bx]© + Dr (26. b) 


As mentioned earlier, the vector r is known from the 
almanac, therefore r is known too. Consequently we can 
treat Dr as a deterministic forcing function. In our case r is 
the vector from the spacecraft to the Sun. For the MAP 
trajectory this vector changes very slowly. Therefore we can 
neglect the term in r and use 


b = [bx]co 


(26.c) 


Augmenting Eq. (26.c) with the dynamics of Filter I yields 

5 


Measurement Equations 

The measurement equations of Filter II are quite simple. 
The vector measurement of Eq. (21) is replaced by the 
simple equation 

M 


b m =[0 0 I,] 


+ V. 


(29) 


Eq.(24) stays the same; therefore, the augmented 
measurement equation for Filter II becomes 







CD 



'K 


' 0 

0 

L 

q 

+ 

’ v b" 

CD 


H 

0 

0 



v„ 

|_ m J 


01 



b 


L 03 


OBSERVABILITY TESTS 

Before carrying on an observability test it is necessary 
to specify the kind of filter that was used in this analysis. 
The filter that was used here is the PSELIKA (Pseudo- 
Linear Kalman) Filter. 6 The idea on which this filter is 
based is as follows. Suppose that the non-linear dynamics 
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and measurement equations of a non-linear system can be 
decomposed in the following way 

x = F(x)x + Bu + w (3 1 .a) 

y = H(x)x + v (3 Lb) 

If F(x) is close enough to F(x)and H(x) is close 
enough to H(x)(note that we do not assume that X is 
necessarily close to x) then instead of the unknown 
matrices F(x) and H(x) we can use respectively F(x) 
and H(x) which are known along the estimation process 
and apply the linear Kalman filter algorithm to the 
measurements using the models 

x = F(x)x + Bu + w (31 .c) 

y = H(x)x + v (3 1 .d) 

The observability analysis requires F(x) andH(x).Thex 
for F(x) and H(x) was the simulation truth model state. The 
logic for this choice of x is as follows. If when using F(x) 
and H(x) in the analysis the latter indicates that the 

system is not observable under these favorable conditions 
then there is no sense in trying to estimate the state vector 
X of the system. If, on the other hand, the analysis 
indicates that the system is observable, then although 
observability and therefore estimability are not assured, it 
makes sense to design an estimator and attempt to estimate 
X . The observability of the system models used by the two 
filters was examined using the following analysis. 

Consider O k , the transition matrix, which corresponds to 
F, the dynamics matrix, of either Eq. (10) or Eq. (27). This 
matrix transforms the system state vector at time t k , X k , 

to X k+1 , the state at time t k+1 . If at a certain time-point, 
t m , the initial state, denoted by x 0 , can be computed, 

then, for our purposes, the system is observable. Adopting 
the common approach to proving complete observability of 
a discrete linear system, we express the first m 
measurements as follows 


y j = H j (x J > 


1=0 


x 0 j = 0,l m 


(32) 


where O 0 (x 0 ) = I , then the last equation can be written 
as 


Y = Wx 0 (33. a) 

where 

Y = [yi y 2 • • y m f (33. b) 

and 


H 0 (x 0 ) 

H,(x,)<t) 0 (x 0 ) 

H 2 (x,)<X>,(x,)<J> 0 (x 0 ) 


(33.c) 


L H n 1 ( X J < 1 > n,- 1 (X m .,)--O 0 (X 0 )j 

Note that in Filter II, H is not a function of the state (see 
equations 23 and 30). That is, the measurement equation 
is linear from the start. If there are n independent rows in 
the right hand side matrix (the observability matrix) in Eq. 

(32) where n is also the number of states, then X 0 can be 
solved for and the state is observable. The observability 


Filter] 



Fig. 4: Rank Development of the Observability Matrix of 
Filter I. 

test of Filter I reveals that the rank of the observability 
matrix for this filter becomes 7 (see Fig. 4), and since the 
filter state is 7 too, the system is observable. It means that 
there is a good chance that Filter I can estimate both the 


Filter II 



Fig. 5: Rank Development of the Observability Matrix of 
Filter II. 
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missing angular rate component as well as the attitude of the 
spacecraft when one gyro fails. We note though that the 
measurement equation of Filter I is only pseudo-linear and 
not truly linear. This fact may cause a problem even though 
the observability test indicates that the state is completely 
observable. The observability matrix of Filter II becomes 6 
and stay at 6 (see Fig. 5) whereas the size of the state vector 
is 9. This disparity indicates that the system is not 
completely observable. It means that we cannot fully 
estimate the missing rate component and the attitude. 
However, maybe their estimated values are close enough to 
the truth. Therefore both filters had to be tested. Simulation 
results of the two filters are presented in the next section. 

SIMULATION TEST RESULTS 

Simulations were run for different initial conditions and 
noise levels. Plots that resulted from a typical run of Filter I 
are presented in Figs. 6-11. For this run the initial y angle 



Time [sec] 

Fig. 6: The x-component of the angular velocity, its 
estimate, and its estimation error for Filter I. 



Time [sec] 


Fig. 8: The z-component of the angular velocity, its 
estimate, and its estimation error for Filter I. 

was 60 degrees, 5 was 35 degrees and r\ was 0. The 
initial estimates were zero for all three angles. The true 
Euler-angle rates were 2-(2 ti)/ 3600 rad/sec for y, 

zero for 5 and 0.928 • (27t)/60 rad/sec for f] . 

According to Eq. (1), using these initial rates and angles, 
the true initial angular rate components were as follows: 

C0 V =0.002 rad/sec C0 V = 0 rad/sec and C0 7 =0.094 

rad/sec. The initial estimated values of the angular 
velocity were zero for all three components. The one- 
sigma gyro measurement noise was 0.01 deg/sec and that 
of the sun sensor was 0.25 degrees. In each figure we see 
the true state superimposed on the estimate of the state. 
Underneath this plot there is a plot of the estimation error. 
In this run y is estimated quite well, but this is not always 
the case. The observability of this angle is quite poor. 




Time [sec] 

Fig. 7: The ^-component of the angular velocity, its 
estimate, and its estimation error for Filter I. 



Time [secj 

Fig. 9: The y angle, its estimate, and the estimation error 
for Filter I. 
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Fig. 10: The 8 angle, its estimate, and the estimation error 
for Filter I. 


- i 

1 True 

1 — Estimated 

-i . i . . 

* : 1 1 : ■ - - 

0 20 40 60 

80 100 120 140 160 ISO 2C 

j 

. \ r i i i 

1 ] i 

; I j | - 


. 


- 50 1 1 1 ‘ 1 1 1 1 1 1 ' 

0 20 40 60 SO 100 120 140 160 180 200 

Time [see] 


Fig. 12: The y angle, its estimate, and the estimation 
error for Filter II. 




Fig. 11: The r| angle, its estimate, and the estimation error 
for Filter I. 
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Fig. 13: The 8 angle, its estimate, and the estimation 
error for Filter II. 


Filter II was also run but, as expected from the observability 
test, the three attitude angles were not observable. The 
estimates of the angular rate components resemble those 
obtained using Filter I but the estimates of the three angles 
were quite poor as seen in Figs. 12-14. 

CONCLUSIONS 

In this paper we analyzed two algorithms for estimating the 
attitude and rate of the MAP spacecraft after one of its two 
two-degree of freedom and one of its vector-measuring 
devices fail. Both algorithms implement Kalman filters for 
estimating its attitude and rate. A special observability 
analysis that was carried out revealed that one of the filters, 
Filter I, was completely observable whereas the other filter, 
Filter II, was not completely observable. It was confirmed 
through simulations, which used the MAP spacecraft 
attitude and rate profile at the L2 Lagrange point, that Filter 
II was indeed unable to estimate the attitude. However, it 




Fig. 14: The T| angle, its estimate, and the estimation 
error for Filter II. 

8 


American Institute of Aeronautics and Astronautics 









was also found that, contrary to the results of the 
observability analysis, when using Filter I the estimation 
quality of one of the attitude angles depended on the initial 
errors. A possible explanation to this contradiction is based 
on the fact that the observability matrix that was used to 
examine the observability of the filter was state dependent. 
Therefore only when we ran the test using the correct state; 
namely, the nominal angular rate vector and Euler angles, 
was the system found to be observable. In reality, however, 
the existing state vector was the estimated state vector. 
Therefore as long as the initial estimate was close enough to 
the true state, the simulation conditions matched those of the 
analysis and the observability test results predicted well the 
performance of the filter. On the other hand, when the initial 
estimate of the state vector was far from the true vector the 
analysis was not valid anymore. The main conclusion of this 
analysis is that when only one vector-measuring sensor is 
active, MAP dynamics equations cannot assist in turning the 
system into a completely observable system. Therefore 
complete attitude estimation is impossible under the 
conditions described in the paper. 
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